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Abstract 

Background: Oxyspirura petrowi appears to be emerging as a nematode parasite that could negatively impact 
Northern Bobwhite quail individuals and populations within Texas and other regions of the United States. Despite 
this eye worm's potential importance in the conservation of wild quail, little is known about the general biology 
and genome composition of 0. petrowi. To fill the knowledge gap, we performed a small scale random genome 
sequence survey, sequenced its 18S rRNA and the intergenic region between the 18S and 28S rRNA genes, studied 
its phylogenetic affinity, and developed a PCR protocol for the detection of this eye worm. 

Results: We have generated -240 kb of genome sequence data derived from 348 clones by a random genome 
survey of an 0. petrowi genomic library. The eye worm genome is AT-rich (i.e., 62.2% AT-content), and contains a 
high number of microsatellite sequences. The discovered genes encode a wide-range of proteins including 
hypothetical proteins, enzymes, nematode-specific proteins. Phylogenetic analysis based on 18S rRNA sequences 
indicate that the Spiruroidea is paraphyletic, in which Oxyspirura and its closely related species are sisters to the 
filarial nematodes. We have also developed a PCR protocol based on the ITS2 sequence that allows sensitive and 
specific detection of eye worm DNA in feces. Using this newly developed protocol, we have determined that -28% 
to 33% of the fecal samples collected from Northern Bobwhites and Scaled Quail in Texas in the spring of 2013 are 
0. petrowi positive. 

Conclusions: The 0. petrowi genome is rich in microsatellite sequences that may be used in future genotyping and 
molecular fingerprinting analysis. This eye worm is evolutionarily close to the filarial nematodes, implying that therapeutic 
strategies for filariasis such as Loa loa would be referential in developing treatments for the Thelazoidea parasites. Our 
qPCR-based survey has confirmed that 0. petrowi infection is of potential concern to quail managers in Texas. 



Background 

Oxyspirura petrowi is a spirurian nematode (Order 
Spirurida) that infects the eyes of quail and other birds [1]. 
In Texas, a 47-56% prevalence has been reported in 
Northern Bobwhites (Colinus virginianus) and Scaled 
Quail {Callipepla squamata) [2-4], Similar infections 
caused by this genus of parasites have also been reported 
in other animals including poultry and zoo animals, where 
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some of them were described as ocular oxyspiruriasis or 
oxyspirurosis [5-10]. Given that bobwhites are experien- 
cing long-term declines throughout their range in North 
America, there is a recognition that populations are 
declining even where suitable habitat conditions exist (e.g., 
Rolling Plains ecoregion of Texas), thereby raising con- 
cerns that parasites such as O. petrowi may be a contribut- 
ing factor (e.g., see a more detailed description at http:// 
www.quailresearch.org). It is likely that infection may 
cause host eye damage and physically impair vision, mak- 
ing birds less competitive in feeding and more susceptible 
to predators (Figure 1). 

Although the eye worm has been considered as a pos- 
sible contributing factor for the decline of wild quail 
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Figure 1 Oxyspirura petrowi adult worms in the eye of a 
Northern Bobwhite collected in Texas in February, 2013 
demonstrating their potential to cause visual obstruction in 
addition to a pathological response resulting from infection. 



populations in the Rolling Plains, little is known of the 
parasite s biology, particularly at the molecular and gen- 
omic levels (i.e., no molecular data were available in the 
GenBank databases prior to this study). Previous know- 
ledge on the relationship of this parasite with other 
nematodes was solely acquired by morphology, which 
also needs to be validated at the molecular level. In fact, 
only a single nucleotide sequence is present in the da- 
tabase for the whole genus Oxyspirura (i.e., a 689-bp 
partial rRNA gene from O. conjuctivalis [GenBank: 
EF417873]). The lack of molecular data severely ham- 
pers our efforts in studying molecular epidemiology and 
transmission routes of O. petrowi, which may be useful 
for developing effective strategies to treat and control 
ocular oxyspiruriasis in wild quail. 

To fill the knowledge gap, we have performed a small- 
scale genome sequence survey (GSS) that provides the 
first batch of genomic sequence data for this nematode. 
Additionally, we have cloned the 18S rRNA, internal tran- 
scribed spacer 1 (ITS1), 5.8S rRNA, ITS2 and partial 28S 
rRNA genes. The small random GSS effort rapidly gener- 
ated -240 kb of sequence information that provided not 
only a snapshot of the quail eye worm genome, but also a 
large amount of microsatellite sequences for future geno- 
typing and population genetic analysis. Based on the newly 
obtained rRNA sequences, we determined the evolution- 
ary affinity of Oxyspirura with other nematodes among 
spirurians with available 18S rRNA sequences. We also 
designed specific primers based on the ITS2 sequences, 
and performed real-time quantitative PCR (qPCR) -based 
molecular detection of O. petrowi from DNA extracted 
from fecal samples collected from Northern Bobwhite and 
Scaled Quail in Texas. 

Understanding the biology of this parasitic nematode 
at molecular levels will enable us to effectively determine 



the prevalence by detecting parasite-specific DNA in 
feces, as well as to identify infected intermediate hosts 
that is otherwise difficult (if not impossible) based on 
morphology of larval stages. Molecular tools would 
enable further study of potential drug targets and target- 
based drug discovery to treat this important nematode. 

Methods 

Isolation of genomic DNA and genome sequence survey 

Adult O. petrowi worms were isolated from the eyes of 
Northern Bobwhites collected in Texas as part of a 
3-year integrated research project called Operation Idio- 
pathic Decline, which was initiated to further our under- 
standing of potentially pathogenic parasites occurring 
within the Rolling Plains Ecoregion of Texas and 
western Oklahoma. All animal experiments were 
performed in accordance with procedures approved by 
the Institutional Animal Care and Use Committee of 
Texas A&M University (protocol # 2011-193). 

After microscopic examination for species validation, 
four worms were rinsed with PBS, placed in lysis buffer of 
the DNeasy Blood & Tissue Kit (Qiagen Inc., Valencia, 
CA), and grinded with a plastic microtube grinder. 
Genomic DNA (gDNA) was isolated from the ground 
worms according to manufacturers protocol for animal 
tissues. For the construction of a genomic library, gDNA 
was first subjected to whole genome amplification using a 
GenomePlex Complete Whole Genome Amplification 
(WGA) kit according the manufacturers standard proto- 
col (Sigma-Aldrich Co., St. Louis, MO). Amplified gDNA 
products were fractionated in agarose gels and fractions 
containing fragments between 0.5 - 2 kb were collected 
and purified using a Gel Extraction Kit (Omega Bio-Tek, 
Norcross, GA). After an incubation at 72°C for 20 min in 
a regular PCR reaction buffer to add a single adenine over- 
hang to the 3'-end, the products were ligated into pCR2.1- 
TOPO vector using a TOPO-TA cloning kit (Invitrogen, 
Life Technologies, Grand Island, NY). 

After transformation, Escherichia coli OneShot TOP10F 
chemically competent cells (Invitrogen) were plated onto LB 
plates spread with 40 uL of 40 mg/mL X-gal and 5 uL of 
200 mM/mL IPTG, and incubated at 37°C overnight. Bac- 
teria from a single white colony were collected into 10 uL 
Milli-Q water in a microtube, from which 2 uL of suspen- 
sion was used directly as template in PCR reactions to deter- 
mine the presence of inserts using a pair of Ml 3 forward 
and Ml 3 reverse primers. Colonies containing inserts with 
desired sizes were further incubated in LB broth containing 
50 ug/mL kanamycin. Plasmids were prepared using a Plas- 
mid DNA Mini kit (Omega Bio-Tek), and sequenced by the 
Sanger dye-terminator sequencing method at the Gene 
Technology Laboratory at Texas A&M University. 

Vector and adaptor sequences were removed using a 
cross-match algorithm, and long inserts were assembled 
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using the Phrap method implemented in the MacVector 
program (version 12.7 A) (http://www.macvector.com). 
All sequences were used as queries to search the non- 
redundant protein and nucleotide databases at the National 
Center for Biotechnology Information (NCBI) by the BLA 
STN, BLASTX and TBLASTX algorithms using the Kori 
Blast program (version 3.4) (http://www.korilog.com). Add- 
itional annotations were performed using the Blast2GO 
program (http://www.blast2go.com/b2ghome), which in- 
cluded InterProScan for identifying protein domains and 
gene ontology (GO) analysis. GO_slim was performed at 
the CateGOrizer server (http://www.animalgenome.org/ 
cgi-bin/util/gotreei) [11]. Contigs were also mapped onto 
the metabolic pathways at the Kyoto Encyclopedia of Genes 
and Genomes (KEGG) using the KEGG Automatic Anno- 
tation Server (KAAS) (http://www.genome.jp/tools/kaas/) 
[12]. Candidate tRNA sequences were examined at the 
tRNAscan-SE server (http://lowelab.ucsc.edu/tRNAscan- 
SE/) [13]. Microsatellite sequences (also known as simple 
sequence repeats, SSRs) were identified using the Phobos 
(version 3.3.12) program (http://www.ruhr-uni-bochum.de/ 
spezzoo/cm/cm_phobos.htm), in which only perfect mat- 
ches with a minimal length of 8 nt and a minimal score at 
8 were reported. 

Molecular cloning of parasite ribosomal RNA 
(rRNA) genes 

The 18S rRNA gene and downstream ITS1, 5.8S rRNA 
and ITS2 regions from O. petrowi were cloned by PCR 
using two pairs of primers: 1) nemal8S_F01 (5'-CCA 
TGC AWG TCT AWG TTC AAA-3') and nemal8S_R01 
(5'-GGA AAC CTT GTT ACG ACT TTT G-3') for the 
nearly whole 18S region; and 2) nemal8S_F1400 (5'-GTC 
TGT GAT GCC CTT AGA TG-3') and nema28S_R68 
(5'-TTA GTT TCT TTT CCT CCG CTT A-3') for the 
region between the 18S and 28S rRNA genes. PCR was 
performed using a JumpStart REDTaq ReadyMix PCR 
Reaction kit containing hot-start high-fidelity DNA poly- 
merase ( Sigma- Aldrich). After treating with regular Taq 
DNA polymerase at 72°C for 10 min, PCR amplicons were 
similarly cloned into the pCR2.1-TOPO vector as de- 
scribed above. At least 10 independent clones from each 
reaction were sequenced, and all reads were assembled by 
Phrap as described above. Regions representing 18S, ITS1, 
5.8S, ITS2 and partial 28S sequences were determined by 
Rfam (http://rfam.janelia.org) [14]. 

Phylogenetic reconstructions 

The assembled O. petrowi 18S rRNA sequence was used 
as a query to search and identity nematode orthologs 
from the NCBI nucleotide databases. Up to 1,000 gene 
sequences were initially retrieved, subjected to multiple 
sequence alignments using the MUSCLE program (ver- 
sion 3.8.31) (http://drive5.com/muscle/) [15,16], followed 



by simple phylogenetic reconstructions by the neighbor- 
joining (NJ) method using the Tamura-Nei nucleotide 
substitution model. Based on the alignment and NJ 
trees, short or identical sequences were individually re- 
moved, and the same procedure was repeated until a 
balanced dataset containing 111 sequences representing 
all major nematode taxonomic groups were identified. 
The dataset was subjected to phylogenetic reconstruc- 
tions by Bayesian inference (BI) using MrBayes (version 
3.2) (http://mrbayes.sourceforge.net) and the maximum 
likelihood (ML) method using TreeFinder (version 2008) 
(http://www.treefinder.de) [17,18]. This approach deter- 
mined the phylogenetic relationships among major taxo- 
nomic groups, in which O. petrowi was placed within 
the spirurians, but the relationship among spirurians 
was not well resolved. Therefore, we resampled the se- 
quences to include only taxa within Spirurida and 
Ascaridida as these two groups displayed a sister rela- 
tionship by this study and previous analyses [19,20]. This 
also allowed us to include more taxa within these two 
groups. The second dataset contained 112 taxa with 
1,544 nucleotide positions and was subjected to phylo- 
genetic reconstructions using BI and ML methods. To 
further resolve the O. petrowi position, we also compiled 
a third dataset containing only taxa with close relation- 
ship with O. petrowi. This small dataset included only 35 
taxa with 1,599 nucleotide positions, and was also 
subjected to BI and ML analyses. 

In all datasets, gaps were removed and only positions 
that could be unambiguously aligned were used in subse- 
quent phylogenetic analyses. In the BI analysis, 1.5 million 
generations of searches for the first and second datasets 
(or 1.0 million generations for the smaller third dataset) 
were performed with 4 independent chains running. 
Searches reached convergence as determined by the aver- 
age standard deviation (SD) of split frequencies reaching 
< 0.01, and potential scale reduction factor (PSRF) values 
for various approaching 1.0 [21]. Bootstrapping ML ana- 
lyses were derived from 200 replicated sequences. In both 
BI and ML methods, the general time reversal (GTR) nu- 
cleotide substitution model was used with the consider- 
ation of fraction of invariance and 4-rate of discrete 
gamma (i.e., GTR + F inv + r 4 ). Majority rule consensus 
trees were visualized using FigTree (version 1.4), followed 
by tree annotations using Adobe Illustrator CS4. 

Molecular detection of 0. petrowi 

Sequence comparison of the rRNA regions between O. 
petrowi and other nematodes indicated that 18S rRNA se- 
quences were less suitable for designing species -specific 
primers, as they were highly conserved among nematodes. 
We hence designed primers based on the ITS2 region se- 
quences for specific molecular detection for O. petrowi: 
QEW_2417F (5'-GGA TTT GCA AGA ATT GTT TCC-3') 
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and QEW_2578R (5'-AAC GTT ATT GTT GCC ATA 
TGC-3') with a predicted product size of 162 bp. To ensure 
desired specificity of detection, we used DNA samples iso- 
lated from Aulonocephalus pennula, a cecal worm collected 
from bobwhite quail in northern Texas as negative control 

We applied real-time quantitative PCR (qPCR) to de- 
tect eye worm DNA from fecal samples from Northern 
Bobwhite and Scaled Quail in Texas. Feces from individ- 
ual or pooled birds were collected at Rolling Plains Quail 
Research Ranch (RPQRR) in Fisher County, Texas in the 
Spring, 2013 via the seasonal trap-and-release program 
in a separate conservation research project. Feces were 
mixed, weighed and placed in lysis buffer included in 
the QIAamp DNA Stool Mini Kit (Qiagen). After one 
freeze/thaw cycle in liquid nitrogen, samples were ho- 
mogenized with glass beads in a Mini-Beadbeater-16 
(BioSpec Products, Inc., Bartlesville, OK) at full-speed 
for 4 min, followed by DNA isolation according to the 
manufacturers protocol for the stool DNA isolation kit. 
A SYBR-green-based qPCR detection was performed in 
20 uL reactions containing iQ SYBR Green Supermix 
reagent (Bio-Rad, Hercules, CA), the QEW_2417F/ 
QEW_2578R primer pair (each at 100 nM) and 2 uL 
stool DNA in a Bio-Rad iCycle iQ Real-Time PCR 
Detection System. The reactions started with sample 
denaturation at 95°C for 5 min, followed by 45 thermal 
cycles at 95°C for 30 sec, 56°C for 30 sec, and 72°C for 
30 sec. Specificity of detection was confirmed by melting 
curve analysis and agarose gel electrophoresis of selected 
PCR products. Selected individual or pooled PCR prod- 
ucts were also sequenced to confirm their identities. 

Nucleotide sequence accession numbers 

Nucleotide sequences generated in this study were de- 
posited into the GenBank database with accession num- 
bers [GenBanlcKFl 10799] and [GenBanlcKFl 10800] for 
rRNA sequences containing type 1 and type 2 ITS1, and 
[GenBank:KG007611] to [GenBank:KG007945] for GSS 
sequences. 

Results and discussion 

Characterization of the O. petrowi genome 

In this small genome sequence survey, we have 
obtained valid sequences for 354 clones. Among them, 
six sequences were determined to be bacterial contami- 
nants, representing 1.6% of the sequenced clones. 
There were no contaminants from birds and other or- 
ganisms, suggesting that the prepared O. petrowi gen- 
omic libraries were of high quality. The limited 
bacterial sequences (top hits were mainly Ralstonia and 
Caulobacter vibrioides) were likely derived from com- 
mensal bacteria in O. petrowi. 

The remaining 348 valid O. petrowi sequences re- 
sulted in 237,239 bp of genomic sequences, which 



ranged from 81 to 1,220 bp (median length = 706 bp) 
for individual contigs. The scale of the survey was rela- 
tively small, but it was sufficient to provide a first snap- 
shot of the genome features for this nematode. The O. 
petrowi genome was generally AT-rich with GC con- 
tents ranging from 18% to 64% for individual contigs 
(overall mean of GC content = 37.8%, vs. 56% to 70% 
for the six bacterial contaminants). Various BLAST 
searches yielded no hits for 137 contigs (-39%), sug- 
gesting that these sequences might be unique to 
Oxyspirura or closely related species. Among the 211 
contigs with hits from the databases, 193 (91.5%) had top 
hits on nematode sequences, particularly on those from 
Loa, Brugia and Wuchereria, but only very few on Ascaris 
or Caenorhabditis, which was congruent with the evolu- 
tionary relationship of Oxyspirura with filarioidea, As- 
caridida and Rhabditida (see Additional file 1: Table SI for 
a complete list of all contigs with annotations and corre- 
sponding BLAST top hits). 

By combining BLAST with InterProScan searches, more 
than half of the contigs with hits were able to be assigned 
into major functional categories (i.e., 121 out of 211 
contigs) (Figure 2). The functionally undefined 90 contigs 
were mainly hypothetical proteins with some containing 
low complexity sequences. Among the 121 annotatable 
contigs, the largest group was enzymes (total 40) that were 
separated into general enzymes involved in various meta- 
bolic pathways (n = 30) and those involved in protein 
metabolism such as protein kinases (n = 8) and proteases 
(n = 2). Examples of enzymes included glycogen synthase 
(contig QEW_195), glycosyl-transferase (QEW_224), his- 
tone acetyltranferase (QEW_156), and succinate dehydro- 
genase (QEW_315); while those for protein metabolism 
included a Ulpl protease family member (QEW_129), as 
well as protein kinases in the AGC/NDR (QEW_74), 
CAM/CAMKL/NUAK (QEW_372) and CK1/WORM6 
(QEW_249) families. Fourteen contigs encoded proteins 
involved in DNA/RNA metabolism, e.g., splicing factor 1 
(QEW_306), DNA replication licensing factor mcm-6 
(QEW_340), and protein containing double-stranded 
RNA binding motif (QEW_379). There were 12 contigs 
containing genes encoding extracellular membrane (ECM) 
proteins, in which 9 contigs were associated with the 
nematode-specific cuticle formation, such as cuticle colla- 
gen precursors (QEW_58, QEW_59 and QEW_135), 
cuticulin-1 (QEW_104), and nematode cuticle collage do- 
main containing proteins (e.g., QEW_80 and QEW_386). 
Other groups include those involved in ribosomal biogen- 
esis (n = 11), molecular interactions (n = 8), ion trans- 
porters (n = 6), cytoskeletal proteins (n = 6), membrane 
proteins associated with cell adhesion (n = 3), and those 
involved in gene expression (n = 2). 

Interestingly, there were three contigs encoding ne- 
matode-specific major sperm protein (MSP), which was 
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Figure 2 Classification of Oxyspirura petrowi genes discovered by the random genome sequence survey by major functional groups 
(A) or gene ontology (GO) terms (B). A list of gene contigs with annotations is provided in Additional file 1: Table S1. 



* NNNA = 

nucleobase, nucleoside, 
nucleotide and nucleic acid 

** metabolites/energy = 
generation of precursor 
metabolites and energy 



grouped together with a fic protein under the "cell de- 
velopment" category in Figure 2. There were also three 
contigs derived from the mitochondrial genomes, in- 
cluding gene fragments encoding cytochrome b, cyto- 
chrome c oxidase subunit IV (COX-IV) and NAHD 
dehydrogenase subunit 5. Finally, 4 contigs were found 
to contain retroelements, such as tigger transposable 
element-derived protein 1-like proteins (QEW_112 and 
QEW_119) and retrotransposon protein (QEW_172). 

Microsatellite sequences identified by the genome survey 

One major purpose of the study was to identify micro- 
satellite sequences (also known as simple sequence 
repeats, SSRs) for future genotyping and population gen- 
etics analyses on this quail eye worm. In fact, O. petrowi 
appears to be rich in microsatellites, in which a total of 
335 units of perfect SSRs were identified with a minimal 
length of 8 nt (Table 1). These included mononucleo- 
tides (228 units), dinucleotides (30), trinucleotides (56), 
tetranucleotides (11) and 10 repeats with 5-8 nucleo- 
tides. At least 98 contigs contained two or more SSRs, 
and 31 contigs contained 3-6 SSRs (Table 1). Examples 
included QEW_123 with 5 for mono-, tri- or tetra- 
nucleotide SSRs; QEW_126 with 5 mono-, tetra- or 



octa-nucleotide SSRs, and QEW_203 with 6 di-, tri- or 
penta-nucleotide SSRs (see Additional file 2: Table S2 
for a complete list of detected microsatellite sequences). 
We also looked at the distribution of microsatellites with 
repeat units of >2 nt, which revealed ~2 or -1.5 times 
more microsatellite sequences are present in contigs 
with no hits in BLAST/InterProScan searches (19.0%) or 
with hits but unknown function (14.4%) than in the 

Table 1 Statistics on the lengths of repeat units and 
numbers of microsatellite sequences per contig in 
Oxyspirura petrowi identified by the genome sequence 
survey 



Length of 
repeat unit 


Counts 


No. microsatellites 
per contig 


Counts 


1 


228 


1 


86 


2 


30 


2 


67 


3 


56 


3 


17 


4 


11 


4 


7 


5 


2 


5 


6 


6 


6 


6 


1 


8 


2 


>7 


0 


Total microsatellites 


335 


Average no. per contig 


1.82 
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annotatable contigs (9.9%) (Table 2). In summary, the 
eye worm genome contains a rich number of microsatel- 
lite sequences with the potential to be further validated 
as potential genetic markers. 

Phylogenetic position of O. petrowi based on 1 8S 
rRNA genes 

Our first phylogenetic analyses based on a large 18S 
rRNA dataset with BI and ML methods produced trees 
that agreed with those produced by others. While O. 
petrowi was clustered within the Spirurida clade, it was 
close to a branch consisting of Tetrameres fissipina and 
an unknown Onchoceridae species. This was likely a re- 
sult caused by a long branch attraction (LBA) artifact 
based on the unusual long branch formed by T. fissipina 
and the Onchoceridae species, as well as by the obvious 
high numbers of nucleotide substitutions in these two 
sequences (data not shown). We also observed potential 
sequencing mistakes for the long 18S rRNA sequence of 
Thelazia lacrymalis (DQ503458). Therefore, we re- 
moved these three sequences from subsequent analyses. 
Because Spirurida was consistently placed as a sister 
group with Ascaridida, as observed in this study and by 
others [19,20], our second dataset contained only se- 
quences from these two groups, which allowed more re- 
liable alignment with more nucleotide positions to better 
resolve the position of O. petrowi within Spirurida using 
Ascaridida as outgroup. Gnathastoma sequences were 
also excluded from the second dataset, as they have been 
shown to be seperate from the rest of the spirurids 
[19,20]. 

Both BI and ML trees inferred from the second dataset 
distinctly separated Ascaridida from Spirurida (Figure 3A). 
Within the Spirurida clade, Dracunculoidea and Camal- 
lanoidea formed two major sister branches, whereas the 
third branch comprised of the remaining families including 
Spiruroidea, Acuarioidea, Physalopteroidea, Filarioidea, 
Habronematoidea and Thelazioidea. Further phylogenetic 
analysis based only on sequences from the third branch 
produced similar tree topology, but with slightly better 
resolution and statistical support (Figure 3B). Acuarioidea, 
Physalopteroidea, Filarioidea and Habronematoidea were 
monophyletic, whereas Spiruroidea was paraphyletic, 



Table 2 Number of microsatellites (SSR) with unit 



length >2 by functional groups* 



Group 


No. contigs 


SSR (unit > =2) 


Percentage 


Annotatable 


121 


12 


9.9% 


Function unknown 


90 


13 


14.4% 


No hits 


137 


26 


19.0% 


Total 


348 


51 


14.7% 



* See Additional file 2: Table S2 for a complete list of microsatellite sequences. 



intermixed with other families. Among them, O. petrowi 
was clustered with Streptopharagus and Spirocerca, which 
in turn formed a sister branch to the Filarioidea, albeit with 
low posterior probability and bootstrap proportion support 
(Figure 3B). At the moment, more sophisticated phylogen- 
etic analyses were unachievable due to the lack of more 
sequences from closely related species, and the lack of 
sufficient sequence data such as the mitochondrial ge- 
nomes and proteins within Spirurida, particularly among 
Thelazioidea. Nonetheless, our study revealed that The- 
lazioidea, including quail eye worm, was closely related to 
filarial nematodes, which implies that therapeutic strat- 
egies for filariasis such as those for L. loa might be refer- 
ential in developing treatments for the Thelazoidea eye 
worms. 

Feature of internal transcribed regions and molecular 
detection of O. petrowi 

In addition to the nearly complete 18S rRNA gene, we 
have also determined the complete sequences of the 
ITS1, 5.8S rRNA and ITS2 regions. Two types of ITS1 
sequences have been observed (GenBank accession 
numbers: KF1 10799 for type 1 and KF1 10800 for type 2 
ITS1), in which the type 1 sequence contained a 17- 
nucleotide insertion at positions 1911 - 1927 (Figure 4). 
Additionally, there were two nucleotide differences be- 
tween the two types of ITS1 sequences at positions 2112 
and 2188 (in reference to KF1 10799). The insertion/de- 
letion was further confirmed by PCR using primers 
flanking the region and subsequent sequencing of PCR 
amplicons. 

While the 5.8S rRNA gene shared significant sequence 
homology to those of other nematodes available the 
NCBI databases (data not shown), both the ITS1 and 
ITS2 regions displayed high sequence diversity (i.e., no 
significant hits in BLASTN searches of the NCBI nu- 
cleotide databases). Because of the insertion/deletion at 
the ITS1 region, we initially designed two pairs of 
primers based on the ITS2 sequence for the potential of 
using nested PCR-based molecular detection of O. 
petrowi: an external primer pair QEW_2373F (5'-AAG 
AAT GTA ATG TTG TGG AGC-3') and QEW_2681R 
(5'-GTA ATC ACA TTT GAG TTG AGG-3'), and an 
internal primer pair QEW_2417F and QEW_2578R (as 
described in the Methods section) that would give 
309 bp and 162 bp products, respectively. However, after 
our pilot experiments indicated that nested PCR was un- 
necessary, we used regular qPCR reactions with the 
QEW_2417F and QEW_2578R primers in subsequent 
detection of O. petrowi DNA from wild bobwhite fecal 
specimens collected in Texas in February - March, 
2013. The specificity of the QEW_2417F and 
QEW_2578R primer pair for O. petrowi was also con- 
firmed by its inability to produce products from DNA 
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Figure 3 Phylogenetic relationship of Oxyspirura petrowi within the Spirurida nematodes as determined by Bayesian inference (Bl) and 
maximum likelihood (ML) methods based on 18S rRNA sequences from Spirurida and Ascaridida (112 taxa with 1,544 positions) (A) 
and from species more closely related to Thelazioidea (35 taxa with 1,599 positions) (B). In both approaches, the general time reversal 
(GTR) nucleotide substitution model was used with the consideration of fraction of invariance and 4-rate of discrete gamma (i.e., GTR + F inv + r 4 ). 
Numbers at the nodes indicate posterior probability (Bl) and bootstrap proportion (ML) supporting values. Nodes highlighted by dots were 
supported by >95% in both Bl and ML bootstrapping analyses. Letter "x" indicates nodes supported by <50% in either Bl or ML analysis. 



isolated from the cecal worm A. pennula that is com- 
monly present in wild quail (Figure 5). 

On the other hand, due to the lack of molecular se- 
quences at the ITS regions from very closely related spe- 
cies, we could not firmly conclude that the relatively 
short primers were absolutely specific to O. petrowi. In 
fact, only six ITS1 sequences were present in the 
GenBank database for the superfamily Thelazioidea, in- 
cluding five from Thelazia species and one from O. 
conjuctivalis (accession: EF417873). However, these ITS1 
sequences were highly divergent from each other and 
from those of O. petrowi (i.e., 47.5% - 48.5% identities 
between the O. conjuctivalis and O. petrowi and 26.1% - 



53.2% between the five Thelazia species within the ITS1 
sequences), which indicated that their ITS2 regions were 
likely divergent enough for speciation. To further ensure 
the quality of detection, selected individual or pooled 
PCR products were also sequenced to validate their 
identities. When the qPCR detection system was used, 
melting curve analysis was performed to confirm that 
PCR amplicons showed the same curves as those of 
positive controls. 

Among the 54 fecal DNA specimens, we have detected 
21 (38.9%) positive samples. However, because the speci- 
mens were collected from both individual and pooled 
quail samples derived from 88 birds, direct calculation 



ITSI-type 1: 1888 GTAAATATCTGAAAGGATGAAAAAAGCAACGCTACTTATTACCACACCTGTTTGTTTTGTTGAAAAAAAT 1957 
ITS1-type 2: 1888 GTAAATATCTGAAAGGATGAAAA ACCACACCTGTTTGTTTTGTTGAAAAAAAT 1940 

Figure 4 Comparison of the Oxyspirura petrowi type 1 and type 2 ITS1 sequences at the insertion/deletion site. Their sequences are 
available at GenBank database with accession numbers [GenBank:KF1 10799] and [GenBank:KF1 10800]. 
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Figure 5 Agarose gel (1.5%) electrophoresis illustrating PCR- 
based detection of Oxyspirura petrowi DNA from quail fecal 
samples. Lane M: 100-bp molecular marker; Lanes EW and CW: regular 
PCR using DNA isolated from adult eye worm (EW, 0. petrowi) and 
cecal worm (CW, A. pennulo) isolated from wild quail as positive and 
negative controls (Ctl); Lanes SI - S7: results of real-time qPCR 
detection from selected positive and negative stool DNA samples. 



of positive rate (i.e., 21/54 = 38.9%) was inappropriate as 
pooled positive specimens might contain both positive 
and negative samples. Therefore, we employed two add- 
itional approaches to estimate the prevalence. The first 
approach was to only calculate positive rate from the 39 
samples collected from individual birds (non-pooled), in 
which 13 samples were positive, giving a 33.3% positive 
rate. The second approach employed software written by 
Dr. Brad Biggerstaff as an Excel Add-In (PooledlnfRate, 
version 3.0) (http://www.cdc.gov/westnile/resourcepages/ 
mosqSurvSoft.html), which was originally developed to 
determine positive rates of viral infections in pooled 
mosquito samples using a maximum likelihood estima- 
tion (MLE) algorithm [22]. By applying a bias-corrected 
MLE estimation, we obtained an infection rate of 27.7% 
with lower and upper limits at 18.6% and 38.7%, respect- 
ively (95% confidence interval). Collectively, we conclude 
that -30% (or between 28% - 33%) of the sampled wild 
quail were infected by the eye worm. The actual rate 
might be even higher, as the fecal samples were only col- 
lected once, rather than continuously for several days, and 
the sensitivity of PCR detection might not be maximal due 
to the inhibitory substances commonly present in fecal 
samples as discussed below. 

The detection of O. petrowi DNA in feces allows rapid 
and sensitive detection of the presence of eye worms 
without the need to examine individual birds. However, 
one needs to be aware of the presence of inhibitory sub- 
stances in fecal samples and the difficulties in releasing 
DNA from eggs or encysted larvae. The presence of 
inhibitory substances could be minimized (if not com- 
pletely eliminated) using the tablets included in the 
DNA isolation kits specifically designed for stool sam- 
ples such as the QIAamp DNA Stool Mini Kit (Qiagen). 
Freeze/thaw cycles combined with homogenization with 
glass beads were necessary to break the eggs or encysted 
larvae to ensure the release of DNA. Furthermore, 
nested PCR might also be used by the addition of a 



primary amplification using the external primers QEW_ 
2373F and QEW_2681R to not only improve the sensi- 
tivity of PCR detection, but to also further eliminate the 
presence of inhibitory substances in a second amplifica- 
tion procedure. 

The life cycle of O. petrowi is not well understood, its 
exact intermediate host(s) as well as its migration details 
in quail. The presence of O. petrowi DNA in quail feces 
indicates that eggs laid by female adult worms in the 
eyes could be passed into the digestive system and re- 
leased with the excrement via the naso-lacrimal ducts to 
the gut route, as previously observed in other eye worms 
such as O. mansoni [23]. 

This PCR detection protocol can also be easily 
adapted to identify intermediate host(s) and perform 
surveillance, which is important in developing effective 
strategies to control the transmission of eye worms in 
the field. Our phylogenetic analysis based on 18S rRNA 
sequences indicated that O. petrowi clustered closely 
with Streptopharagus and Spirocerca (Figure 3). It is 
known that birds are a paratenic host for Spirocerca 
lupi infection in their life cycle between dogs and dung 
beetles [24]. Dung beetles are also the intermediate 
host for Streptopharagus [25]. These observations sug- 
gest that dung beetles might be worth examining as 
one of the potential intermediate hosts. Indeed, the de- 
tection of O. petrowi DNA in various insects including 
dung beetles is currently ongoing as part of a separate 
project in determining the intermediate host(s) and 
transmission route(s), and the data will be reported 
upon the completion of the survey. 

Conclusions 

We have performed a small-scale genome sequence sur- 
vey (GSS), which not only rapidly generated a large 
number of molecular sequence data for the first time 
for O. petrowi, but also provided a snapshot of the gen- 
ome for the eye worm in quail. The survey also identi- 
fied a large number of microsatellite sequences that 
may be employed in further genotyping and population 
genetics studies. Our phylogenetic reconstructions 
based on 18S rRNA sequences indicated that Spi- 
ruroidea was paraphyletic, while O. petrowi, Strepto 
pharagus and Spirocerca formed a sister clade to the fil- 
arial nematodes. The obtained ITS sequence data also 
permitted us to design specific primers for molecular 
detection of O. petrowi in fecal samples, which may also 
be adapted to detect this nematode in insect intermedi- 
ate hosts for surveillance and developing strategies to 
control the transmission of eye worms from intermedi- 
ate hosts to quail. We also determined that -28% - 33% 
of the birds were O. petrowi positive, suggesting that eye 
worm was a significant parasite in at least some quail 
ranches in Texas. 
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